WANNIER FUNCTIONS OF ELLIPTIC ONE-GAP 

POTENTIAL 



E D BELOKOLOS, V Z ENOLSKII, AND M SALERNO 



Abstract. Wannier functions of the one dimensional Schrodinger 
equation with cUiptic one gap potential are explicitly constructed. 
Properties of these functions are analytically and numerically in- 
vestigated. In particular we derive an expression for the amplitude 
of the Wannier function in the origin, a power series expansion valid 
in the vicinity of the origin and an asymptotic expansion character- 
ising the decay of the Wannier function at large distances. Using 
these results we construct an approximate analytical expression of 
the Wannier function which is valid in the whole spatial domain 
and is in good agreement with numerical results. 

1. Introduction 

The spectral analysis of Schrodinger operators with periodic poten- 
tials has been investigated since the arbour of quantum mechanics. In 
spite of this, it still represents a non exhausted topic of ever continuing 
interest. From one side, it plays a fundamental role in condensed matter 
physics where it provides the mathematical basis of the quantum the- 
ory of solids. From the other, Schrodinger operators with periodic and 
quasi-periodic potentials play an important role in the integration of 
the Kortweg-de Vries (KdV) equation. Eigenstates of these operators, 
also called Bloch functions (BF), have been extensively studied during 
the past years by several authors (see [111 121 CD] ) ■ An expression of the 
BF in terms of hyperelliptic ^—functions was given in Ref. ^21- These 
studies were further developed in Ref. ^H] where an algebro-geometric 
scheme for constructing solutions of non-linear equations was given in 
terms of the Baker— Akhiezer function. This function is uniquely de- 
fined on the Riemann surface associated with the energy spectrum and 
its properties are natural generalizations of the analytical properties of 
the BF of finite-gap potentials. 

Besides BF, another set of functions which play an equally impor- 
tant role in condensed matter physics are the Wannier functions (WF) 
These functions are related to BF by a unitary transformation 
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and form a complete set of localised orthonormal functions spanning a 
Bloch band. The properties of these functions were first investigated 
by Kohn in 1959 jTSj in a classical paper in which the asymptotic 
decay of Wannier functions was characterised for the case of centro- 
symmetrical one- dimensional potentials. Since then, a large amount of 
work has been devoted to this topic and we mention here results only 
some of them. The projection operator technique was developed for 
construction of the Wannier functions and the Wannier functions were 
studied in the n-dimensional lattices [7j, [S]. The localization problem 
for the Wannier functions was considered in a 1-dimensional case |6j. 
These functions were utilized with success in new practical methods 
for the electron energy calculation of solids (e.g. and in a number 
of modern calculation problems such as the photonic crystal circuits 
|3].The Wannier functions represent the ideal basis for constructing 
effective Hamiltonians of quantum problems involving spatial localiza- 
tions induced by electric and magnetic fields [H],!!!], [TT] . 

In spite of this, the properties of these functions are still not fully un- 
derstood and, except for simplest cases, there are no models for which 
the analytical expression of the WF can be explicitly given. On the 
other hand, the recent results achieved in the field of completely inte- 
grable systems open the possibility to investigate analytical properties 
of the WF. Quite interestingly, WF have not been considered in the 
field of finite gap potentials. 

The present paper represents a first contribution in this direction. 
In particular, we consider WF of Schrodinger operators with one-gap 
potentials and use the well developed theory of elliptic functions to 
investigate their properties with sufficient completeness. As a result we 
derive: i) an exact value for the amplitude of the WF at the localization 
site; ii) an asymptotic expansion characterising the decay of the WF 
at large distances; in) a power series expansion valid in the vicinity of 
the localization site. Using results we construct an approximate 

analytical representation of the WF which is valid in the whole spatial 
domain. These results are shown to be in very good agreement with 
the WF obtained by means of numerical methods. 

The paper is organised as follows. In Section 2 we discuss the basic 
properties of the BF for one gap potentials. In particular we introduce 
basic definitions, discuss the basic properties of BF and derive the 
analytical expression of their normalization constants. Section 3 is 
devoted to the study of the WF. After recalling the basic definitions 
we derive the main results of the paper i.e. points ij-iii) listed above. In 
Section 4 we construct an approximate analytical expression of the WF 
and compare the results of our theory with WF obtained from the basic 
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definition using numerical tools. Finally, in Section 5 we summarise the 
main results of the paper and briefly discuss future developments. 



2. Properties of Block functions of one-gap potential 

2.1. The Schrodinger equation with one-gap potential. In the 

paper we use standard notations and facts of the theory of elliptic 
functions. In particular we use the well known Weierstrass p—, a— 
and (^—functions. 

The periodic elliptic one-gap potential U{x) considered in this paper, 
is expressed in terms of the Weierstrass p— function as 

U{x) = —2p{u), u = ix + UJ, X G M, 
^'^'^^ U{x + nT) =U{x), nGN, 

where T = —2iuj' is the period of the lattice (notice that T is real and 
lX(x) = 2p{ix + uj) is a smooth periodic real function). The Schrodinger 
equation associated with potential (j2.1|) 

(2.2) dl^{x;E) + {E-U{x))'^{x;E) = 0. 

As is well known (see for example [31 and references therein), this equa- 
tion admits eigenf unctions \E'(x;i?) which satisfy the Bloch condition 



(2.3) ^(x-T;E) = e-'^(^)^^(x;E), 
where ii^ is the energy given by 

(2.4) E = p{v), v = a + uj', a G M, 
and k{E) is the quasi- momentum given by 

v' 

(2.5) k{v) = Civ)-—v, v = a + uj', a G M. 

uj' 

Eigenf unctions of type ()2.3|) are called Bloch functions. Notice that 
the dependence of the quasi-momentum on the energy (and vice versa) 
arises from the elimination of the parameter a from Eqs. ()2.4j) . ()2.5|) . 
The BF can be written in explicit form as 

(2.6) <i^{u;v) = C(.;)4^^expK(a)}, 

a[v)a[u) 

or, alternatively, as 

(2.7) ^{u,v) = D{v)y/p{u) - p{v)exp ^ ^'^""^ '^'^ 



p{u) - p{v) 
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where C{v), D{v) are proper normalization constants. In the following 
we shall use both representations for the BF. Notice that the BF con- 
sidered as function of k instead of a, is periodic in the reciprocal space, 
with period 

^~ IT'' 
2uj = —. 

uj' 

BF has the following following periodicity properties 
(2.8) 

^!iu + 2nijj] v) = exp{2nuk(v)}'^(u; v), kiv) = ((v) — v—, 
(2.9) 

m{u + 2n'J] v) = exp{2n'uj'k'{v)}'^{u; v), k'{v) = ({v) - v—^. 

2.2. Normalization of the Bloch function of one-gap poten- 
tials. Since the normalization of the BF plays an important role in 
the construction of the WF (see next section), we shall show how to 
compute the normalization constant, although this question has been 
considered in Chapt. VIII of [2]. 
We normalise the BF according to 

(2.10) 27r(|vl/(x,E)|2) = l, 

where 

1 /-^^ 

(/(a;)) = lim - / /(x)da;. 
The following proposition is valid. 

Proposition 2.1. Normalised Bloch functions of elliptic one-gap po- 
tentials are of the form 

(2.11) 

r]'' 



^(x; a) 



(27r)V 



2 



-1/2 / N 



o[v)a[u, 



exp{i;?7 + (m — uj)C,{v)}, 



where m = ix + cj, x G R; v = a + uo' , a G M. 
Proof. Let us denote the normalised BF as 

where $(u; v) is the non- normalised BF 

a[v)a[u) 
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and C{v) is the normalization constant defined by Eq. ()2.1()j) 



27T 



C{v] 



2u}' 



^(u: v] 



du = 1. 



The complex conjugated (non- normalised) BF is 
a(v - 2uj') ^ ' 



a[v + u) 



-(v-u)')2'q~{u-u))C{v) 



a[v)a[u) 

where the following elementary equalities were used 



a{z) = a{z), C{z) = C{z), {u - uj) = -{u - u), 

u = —u + 2uj, V = V — 2uj' , 

C(t;-2a;')=C(^)-2r/', 

a{v - 2uo') = -a{v) exp{-{v - c^') V), 

cr(M — 2uj) = — cr(n) exp(— (n — uj)2ri). 



By multiplying the above expressions of v) and v) we get 

where in the last step we have used the well known formula 
a{v-u)a{v + u) 

The normalization condition can be then written in the form 

C{vf 



1 = 27r 



2uo' 



27r 



C{v) 



Uu 



-{v-Lu')2r] 



Since 



1 
2^ 



u)+2u)' 



p{u)du = ^[Cico) - Cico + 2cu')] = 

2uj UJ 
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Thus we have obtained for the normahzation constant the expression 



(2.12) 



C{v) 



J9 



(27r)V2 



7] 



UO' 



-1/2 



with an arbitrary phase factor exp(i6'), 6' G M. In the following we fix 
this factor as 

exp(i^) = exp(a;'77 — i(7r/2)). 

□ 

The normalised BF, satisfy a number of useful properties 

under the action of symmetry operations. For centro-symmetrical po- 
tentials the transformation x — > —x of the lattice corresponds to a 
transformation in the Jacobian u ^ u — —u + 2a;, and the following 
propositions can be proved. 

Proposition 2.2. 

^(m; v) = "^{-u + 2uj; v) = ^{u; v). 

Proof. 

^'(m; w) = "^{-u + 2uj;v) 



= i(27r)-^/^ 
= i(27r)-V2 



-Piv) - ^ 



a{v)(7{u — 2uj) 



exp[^;77 -{u- to)C,{v)\ 



/H -1/2 



a{v + u) 
a{v)a{u) 



exp[—vr] — {u — u!)({v)] — '^{u; v). 



□ 



Similarly, that the transformation a —a corresponds to v ^ v 
—V + 2uj', and the following proposition is valid. 

Proposition 2.3. 

^(m; v) = "^{u; -V + 2a;') = '^{u; v). 

Proof. 

*('u;'D) = -v + 2J) 

"^/^ g(-v - m + 2J) 



-i(27r) 



-1/2 



-i(27r)-V2 



a\—v + 2uj')a\u 
X exp[(-i; + 2a;')7^ + (k - a;)C(v - 2a;') 

^^"^ a{y + u) 



a{y)a{u) 



exp[— wr/ — {u — Lo)C,{y)\ 
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□ 

The above propositions can be used to study the elementary proper- 
ties of the BF, k) = \I/(u(a;); v{k)), where u{x) = ix + u and v{k) 
is the inverse of the function k{v) = C(f ) — [r]' /uj')v. To this regard 
note that 

x{u) = —x{u), k{v) = —k{v). 

The following two properties are easily proved. 
Property 1. '^{—x; k) = \l/(a;; k), 

Proof. 

^(-x; k) = ^(m, v) = ^{u; v) = "^{x; k); 

□ 

Property 2. ^(x; -k) = k), 
Proof. 

v[/(a;; -k) = ^(n; v) = "^{u; v) = "^{x; k). 

□ 

3. Analytical properties of the Wannier function of 
elliptic one-gap potentials 

3.1. Definition and basic properties. In 1937 G. Wannier intro- 
duced a complete set of functions for an electron in a lattice structure 
[TBj . The Wannier functions, Wn{x), are defined as 



/ \ 1/2 ^tt/T 

(3.1) ^"(^)=V2^j J ^^^n{x;k)dk. 

where the integral is made on the Brillouin zone. WF for the 
Schrodinger operator with periodic potential U{x), U{x — mT),m G Z 
are localised linear combinations of all the Bloch eigenstates of a given 
n— th spectral band. One can easily prove that if the BF is normalised 
according to Eq. ()2.10|1 . then the WF is normalised on the full line, 

\WJx)\Mx = 1. 



oo 



Using the translation operator one then constructs a countable set of 
WF: Wn\x) := Wn{x — IT), I G Z which is complete and forms an 
ortho-normal basis 

/oo 
W:\x)wiPix)dx = 5nn'5w, I G Z. 
-oo 
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The inverse transformation allows to express a BF in terms of WF as 

/ ^ N 1/2 oo 

(3.2) ^^{x;k)=(-] 5^ iy«(x)e^"^^ 

^ / l=-oo 

In the following we shall omit the band index n since we deal only with 
one band. Properties of WF of one dimensional periodic potentials were 
studied by W. Kohn |13j where he proved that for every band there ex- 
ists one and only one WF which satisfies simultaneously the following 
three properties l)W{x) = W{x)] 2)W{-x) = ±W{x); 3)W{x) = 
0{exp{—h\x\)), where h > 0. In the following we investigate the ana- 
lytical properties of the WF for the one-gap potential in Eq. 1)2.11) . In 
this case the WF is given by the formula 




Using the properties of the BF, "^{x; k), the following basic properties 
of the WF can be proved. 

Proposition 3.1. W{x) = W{x). 

Proof. 

W{x)=i — \ Re J '^{x;k)dk 

= f — J Re y m{x]k)dk = W{x). 

□ 

Proposition 3.2. W{—x) = W{x). 
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Proof. 

ly(-x) = ( — J Re y ^(-x; A;)dA; = f — J Re y ^(x; A;)dfc 
= \ — \ ^^J '^{x;k)dk = W{x). 



□ 

3.2. Power series expansion of the Wannier function at x=0. 

We shall construct in this section the power series expansion of the 
Wannier function of one gap potential. 

Theorem 3.3. The Wannier function of the lower energy band for the 
one gap potential admits the following power series representation 

(3.4) W{x) = -J^^^p^"'^ 

where the coefficients of the expansion ^^.4]) are given by the for- 
mula 



(3.5) W2p = J2^i<lp,i- 

1=0 

Here 

(3.6) M, = ^V^7^t^l^^,-Ue.-e,)'-' 



where F{a,b; c; z) is the standard hypergeometric function, 62,63 are 
branch points of the elliptic curve andqp^i are coefficients of polynomials 
in p{v) 

p 

1=0 

defined by the recurrence 

(3.7) Qp{p{v)) = I 2m - 2 ) <^'»-p-i<5m(p(^)) 

m=0 ^ ^ 

with 



10 



E D BELOKOLOS, V Z ENOLSKII, AND M SALERNO 



First few coefficients of the expansion \3.4\) are 
= Mo, 
W2 = Mi + 2eiMo, 
(3.8) W^ = M2 + 4eiMi + (4e? + 2p"(u;))Mo, 

1^6 = Ms + 6eiM2 + {W{uj) + 12e?)Mi 
+ (2p(^^)(cj) + 28p"(u;)ei + 8e?)Mo. 

Proof. We have 



(3.9) 



1 t/T 

2 



V^^(a;) = ) Re / ^(x; fc)dfc. 



Because the BF \l/(x; A;) is even in x we can write it in the form of 
Taylor expansion 



(3.10) 
where 



, p=l, .... 



If we substitute the Taylor expansion ()3.1()|1 to the ()3.9|) we obtain the 
expansion (|3.4p with the following coefficients 



(3.11) W2p= ' Re y 



Using Schrodinger equation we obtain easily for the ^E'2p(w) a recur- 
rent relation 

^2p(t^) = Y. \ 2^-2 ) <^P-'-i^2/(^) 

0o = 2ei + p(t;), 0p = 2p(^P)(a;). 
The form of this relation leads to conclusion that 

where 



tt/T 



(3.12) 



^'o(^^) = ^(cu;t;) 



1 / p(f ) - ei 
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and Qp{p{v)) arc polynomials of the p-th order in p{v), 

p 

Qp{p{v)) = ^qp,ip\v). 

1=0 

Similarly to ^2p{v) the polynomials Qp{p{v)) satisfy the following re- 
current relation 

(3.13) Qp{p{v)) = '^i 2m -2 ) ^ra-p-lQm{p{v)) 

m=0 ^ ^ 

with 

0o = 2ei + p(^;), 0p = 2p(2^')H. 

In particular the first few polynomials Qp{p{vy) are 

<5i(p(^)) = p(^) + 2ei, 

Q2(p(^)) = + 4eip(^;) + 2p"(a;) + 4e?, 

^^■^^^ Q3(p(^^)) = p{vf + 6eip(T;)2 + (I4p"(a;) + 126?) p{v) 

+ 2p(^^)(o;) + 28eif/H + 8e?. 

Next we calculate the integral expressions of the coefficients W2p- We 
show that the following formula is valid 

Ml = -2 ' J p{vy (^p{v) + ^{u; v)dv 





1 1^' 



/" I rf' 



After the substitution p{v) = s, the computation is reduced to the 
derivation of the complete elliptic integral 



(3.15) M, = /.S/ , 't/ 

^ ^ TT 7 V («-e2)(s-e3) 

By introducing the new variable 

- ^ ~ ^3 

62 - 63 
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the integral Mi acquires the form 



(3.16) = yi^^ y ((62- 63)^-63)' W^^^-^dt, 



where 

(3.17) a;=./^JSZ^, 

is the Jacobi modulus of the elliptic curve 

(3.18) Y^ = {X- e2)(X - 63) + 3 

Using the integral representation of the hypergeometric function (see 
e.g. P) 

Tic) t^-^(i-tY-^-^ 

(3.19) Fia,kc-,.) = 1^ At, 
we obtain the required expression 

(3.20) i='> ' ^ 
xF(-i.; + i;.; + l;p). 

□ 

It is worth to note that the coefficient Wo gives an exact value of the 
amplitude of WF at the localization site x = 0, 



(3.21) Wo = ^^oo'es + v'F (^-^, ^; 1; k'^ . 

This amplitude can also be written in the alternative form 



^21 



(3.22) Wo = — V^'e3 + v'Eik), 

71 

where E{k) is the complete integral of the second kind depending on 
k. Also note that by using the relations 

z-^ + b] = bF{a, b + l;c;z), 



dz 



'1 — z)-^ + c — a — b 
dz 



F{a, b- c- z) = i^-^M^— ^F(a, 6; c + 1; z) 
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one can express all hyperegeometric functions F(— |, j + |; j + 1; k"^) 
in (j3.2Up in terms of the derivatives of F{—^, |; 1; fc^) with respect to 
k"^, and therefore the whole expression can be written in terms of the 
complete integral E{k) and of its derivatives. 

We remark that the quantities p^'^^\uj),j = 1, . . . , can be computed 
in recurrent way, the some first of them are [T] 

3 



p" 




= 3! 








= 5! 


(^■^ 






= 7! 





22 ■ 5^^'^ ~ 2^^' 

1 o 1 



3.3. Asymptotic expansion of the Wannier function. In this sec- 
tion we obtain the asymptotic expression for the WF at x — +oo by 
the steepest descent method. This method (see, e.g. P) permits to 
compute the asymptotic expression of integrals of the type 



F{x) = J f{z)exp{xS{z)}dz, 



where 7 is a contour in the complex plane and the functions f{z) and 
S{z) are holomorphic in the vicinity of 7. When a saddle point zq, 
which is defined by the equation 

^S(.„) = 0. 

does not coincide with the edges of the contour, the asymptotic formula 
of the integral F{x) reads 



2n 



In our case we must use a non-standard variant of the steepest descents 
method, since the exponential in the integrand will have the form xS{z) 
only at — >• +oc and f{zQ) = 0. 

Proposition 3.4. At x ^ 00 the Wannier function of the lower energy 
band for the one gap potential has the following asymptotic expression 
(3.23) 
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where v is a solution of the equation 

(3.24) p(^) = _^, or k'{v)=Q, 

such that the complex number u' k{v) = uj'{({v) — {r]'/uj')v) has negative 
real part. This means that 

W{x) c:^ exp{—h\x\)\x\~^^^, |a;| ^ oo, 

where h = \k{v)\ and v is defined by the equation k'{v) = 0. 

Proof. We use the expression of the WF for the first energy band [es, 62] 
given in Eq. ()3.3|) . Since 

a{v — u) a{u — v) (y{u) aiu — v) 
a{v)a{u) a{ijj — v) a{u) a{v)a{ijj) 

= exp I -CWld.) 

= [p{v)-eiYl^exp\-vr] + j [((s - t;) - C(s)]ds| , 
we have that Eq. ()3.3p can be rewritten as 



(3.25) W{x) = Re £ ^ [e, - piv)]'/' 



dv 



X 



expjy [C(s - t;) - C(s) + C(tO]ds| dv 



When a; ^ +CX3 we can calculate the integral in Eq. ()3.25|) by the 
steepest descents method. To this regard we remark that the argument 
of the exponential, as a function of v, has saddle points which are 
defined by the equation 

''\as-v)-as)+av)]ds = o, 



dv 

or, in other words, by the equation 

. X (ju-v) -({u-v) 
p{v} = . 

U — CO 

At X ^ +00 the last equation attains the form ()3.24p . Since p{v) = 
p{v) we have that if v is saddle point then v is also a saddle point. On 
the other hand, p{v) is an elliptic function of the second order, so it 
takes every value twice in the fundamental domain, this implying that 
there are two saddle points, say f 1, f2, in the fundamental domain. The 
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sum of the values vi,V2 must be a period of the lattice, which in our 
case means that Vi + V2 = 2{uj + u'). It is not difficult to show that 

vi=uj + uj' + il3, V2= oj + uj' -\I3, (3 e M, 

i.e. the two saddle points fi,f2, are situated in the spectral gap. The 
periodicity of the Weierstrass function p{z) in the complex plane give 
rise to countable set V of saddle points, 

V = {vi + 2niLj' , V2 + 277-20;' : rii, n2 G Z}. 

In order to build the proper asymptotic expression for the Wannier 
function W{x) we must select from this set a special saddle point which 
we denote by vq. In the neighbourhood of vo we have 



dk(v) 
dv 



+ 7^{v~v,f [C {s - V,) + C {vo)]ds = / [C(.-t;o)-C(^) + CMd. 



[U — ijJ) — 



u — uo 



{v - v^f 



[C{s - Vo) - C{s) + avo)]ds -{u- uj)-p'{vo){{v - Vo)^. 



Substituting the last two expressions into the integral representation 
of the WF in (jT^ , we obtain 



1/2 



Wo{x) ^ Re <! (^e, + exp { / [C(. - vo) - ({s) + avo)]ds 

X ^ dv[-p'{vo)]'/\v - Vo)'/' exp {-^(« - u;)p'{vo)iv - vo)' 

^ , \/—2iuj' ( vl\'/'^ oivo — u) 

Re < y ei + -7 , \ I \ exp {{u - uj)avo)} 

TT V ^ J cr[vo)(y[u) 

X ^ dr[-p\vo)f\"' exp - uj)p\vo)r' 
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where Cq is a contour passing tlirougli the saddle point vq. The integral 
in the last expression can be calculated as 



Co 



'2p'(vo) 



1/4 



X4 Jo 



e-H-'/^dt 



'2p'{vo) 



r (!) 



3 



Notice that, since p'{vq) = -2[p(vo) -ei]^/2[p('yo) -e2]^/^[p('t;o) -63]^/^ 
and 63 < 62 < pivo) < ei, we have that p'iyo) — — i|p'(i'o)|. 
For the function H^o(2^) we finally obtain 



Wo{x) ~ Re 



\UJ' 



TT 



1/2 



X exp 




2p'(^;o), 

[Q{s - vo) - as) + avoMs 



1/4 



3 

X4 



ei + — 
a;' 



1/2 



=(«-'^)C(fo) 



i 


r (1) 1 


.2p'(i.„). 


Xi J 



(7(t'o)(7(lt) 

The asymptotic behaviour of the function Wq{x) at x — > +00 is defined 
by the factor 

a{vo 



f{u,vo) = 



which satisfies the relation 



u 



a{u) 



exp[{u - uj)Civo)], 



f{u+2nw', Vo) = f(u, Vo) exp[2nw'{C{vo) -vo)] = f{u, vq) exp[2na;'A;(^;o)]. 



The saddle point Vq must be chosen in such a manner that the com- 
plex number u!'k{vo) has negative real part. Prom the previous con- 
siderations it follows that the point vq is also a saddle point. The 
asymptotic behaviour of the function Wo{x) at x — * —00 is de- 
fined by this saddle point which corresponds to the complex number 
cu'kiyo) = Lj'{((vo) — {rj' /uj')vo) with positive real part. □ 

It is appropriate to make here some remarks. 
According to the theorem 



^ exp(-|a;||A;(Q;o)|)|a;| |a;| 



00. 



It is easy to understand such an asymptotic behaviour of the WF at 



\x\ 



00 if we take into account that 



W{x) ^Re{ {k- ko)^e 



\j3Akx 



2sin(/?7r)r(l + /?)x 
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where ko is a branching point of the energy E{k), and that, due to a 
normahzation constant of the wave function, 



^(fc) ~ (k-ko) 



-1/4 



the equahty /3 = — 1/4 is vahd. As far as we know this asymptotic law 
was mentioned for the first time in jB]. 

It is of interest to note also that the equation 

p{v) + ^ = 

UJ' 

has obviously the following solution 

oo 

_l_ . dx 



a/ 4x3 - g2X - gs 

More general problem to solve the equation 

p{v,uj,uj') = c{uj,u'), 

is a well known mathematical problem in the theory of elliptic func- 
tions. Solution of the problem in terms of Eisenstein series is presented 
in the paper 0. 

In the above theorem we have obtained results for the Wannier func- 
tion of lower energy band. Results for the Wannier function of higher 
band, 

(3.26) 



n ./„, V dv (7(v)a(u) 



k{uj) = CM - = ^' 

are similar to ones presented above and as a result of that we omit 
appropriate considerations. 

Let us now discuss two limiting cases. The limit r = u'/u 
corresponds to the free electron case or the empty lattice case. In this 
limit the energy gap is zero, ei = 62, there are no saddle points and as 
a result we have the well known free electron WF 

W{x) = sm — . 

The limit r = u'/u ioo corresponds to the case of tightly bound 
electrons. In this case the width of the lower energy band is zero. 
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-2 -| 1 1 1 1 1 1 1 

-1.4 -0.7 0.0 0.7 1.4 

k 

Figure 1 . Energy bands of the one gap potential with 
parameters of the eUiptic curve ei = 2, 62 = —0.5, 63 = 
-1.5 



62 = 63, the appropriate wave function looks as follows, 

^ /q;\V2 1 ^ 2 

V2/ cosh ax' ° 

where Eq is the binding energy. The wave functions of higher energy 
bands are of a form. 



^(x,k) = ^ , =(|A;| + iatanh(aa:))e^'^^ E = k'^. 

In the next section we shall compare our analytical results with nu- 
merical ones. 
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W(x) 
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0.3 

0.0 



-10 -5 5 10 

X 

Figure 2. The Wannier function associated to the 
lower band of the one gap potential. The branching 
points of the elliptic curve are fixed as in Fig. ^ The am- 
plitude of the function in the origin is W{0) = 0.93. The 
continuous curve denote the exact expression obtained 
from numerical calculations, the dashed line corresponds 
to part of the WF approximated by the asymptotic ex- 
pansion, while the dotted line denotes the part obtained 
from the power expansion near the origin. The arrow 
shows the point where the two different analytical ex- 
pansions are joined. 
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1.2 - 



W(x) 



0.0 



0.8 - 



0.4 - 




-10 



-5 







5 



10 



X 



Figure 3. Same as in Fig|21but for a different set of pa- 
rameters. The branching points of the elhptic curve are 
ei = 6, 62 = —2.0, 63 = —4.0. The value of the function 
in the origin is W{0) = 1.22302. The continuous curve 
denote the exact expression obtained from numerical cal- 
culations, the dashed line corresponds to part of the WF 
approximated by the asymptotic expansion, while the 
dotted line denotes the part obtained from the power 
expansion near the origin. The arrow shows the point 
where the two different analytical expansions are joined. 

4. Approximate analytical expressions of Wannier 
function and numerical results 

The results of the previous section permit to construct the following 
approximate expression of the WF for one-gap potentials: 



W{x) 



Wo + W2X^ + WiX^ + Wqx\ 

for |a;| < Xq 



(4.1) 




for |x| > Xq 
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Figure 4. Asymptotic decay of the Wannier function 
in FigI21 in semi- log scale. The continuous curve repre- 
sents our analytical approximation while the dotted line 
is obtained from direct numerical calculations of the WF. 




Figure 5. Asymptotic decay of the Wannier function 
in Fig. 121 in semi-log scale. The continuous curve repre- 
sents our analytical approximation while the dotted line 
is obtained from direct numerical calculations of the WF. 
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where u = ix + uo ^ v = v+ + uo' ^ the coefficients Wq . . . , Wq are given by 
the formulae (j3.8|) and the point Xq is chosen so to satisfy the normal- 
ization condition = 1. 

To check the validity of this expression we shall compare the WF ob- 
tained from Eq. ()4.1|) with the one obtained directly from the definition 
()3.H) by numerical methods, using the expression of the normalised BF 
in Proposition l2.1l In Fig^we show the band structure obtained for the 
one gap potential with parameter values ei = 2,e2 = —0.5,63 = —1-5, 
while in Fig. |21we depict the WF associated to the lower band. We see 
that the agreement between the analytical approximation and direct 
numerical calculations is excellent both in proximity of the origin and 
far away from it, these being the regions of validity of the correspond- 
ing expansions. In the intermediate region, however, some discrepancy 
appears. It is possible that for some set of potential parameters the 
validity regions of the two expansions (near the origin and far from the 
origin) can overlap at some point Xq so that it is possible to join them 
into the single smooth analytical approximation (j4.1|) which stay close 
to the exact numerically curve in the whole spatial domain. Existence 
of the point of such a kind is possible only for one-gap potential. Such 
a good matching of two expansions is shown in Fig. |21 where the WF 
of the lower band for a potential with another set of parameters is de- 
picted. By comparing Fig. El with Fig. |21we see that the discrepancy 
in the intermediate region is smaller for WF which are more localised. 
This can be understood from the fact that a faster decay of the function 
(see Figs. I4|5I below) allows the asymptotic expansion to work up to 
points which are very close to the origin. In Fig. s|3] and El we show the 
asymptotic decay of the WF depicted in Figs. |2l El respectively from 
which we see that a stronger localization of the function corresponds 
to a faster asymptotic decay. The linear decay observed in the semi-log 
plots of these figures, is fully consistent with the exponential decay of 
the WF of one-gap elliptic potentials predicted by our analysis. 



5. Conclusions 

In this paper we have investigated the properties of the Wannier 
functions of the Schrodinger operator with one gap potentials. As a re- 
sult we have derived the exact value for the amplitude of the functions 
in the origin, as well as, an asymptotic expansion characterising the 
decay of the function at large distances and a power series valid in the 
vicinity of the origin. Using these expansions we have constructed ap- 
proximate analytical expressions of the Wannier functions and showed 
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that they are in good agreement with the ones obtained from numerical 
results. 

We remark that the developed approach can be generalised to the 
case of finite-gap potentials of more complicated type, like elliptic finite- 
gap potentials and general finite-gap potentials. We shall discuss these 
problems in a future publication. 
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